Noise driven dynamic phase transition in a one-dimensional Ising-like model 
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The dynamical evolution of a recently introduced one dimensional model in [I| (henceforth referred 
to as model I), has been made stochastic by introducing a parameter /3 such that /3 — corresponds 
to the Ising model and /3 — > oo to the original model I. The equilibrium behaviour for any value 
of /3 is identical: a homogeneous state. We argue, from the behaviour of the dynamical exponent 
z, that for any /3 7^ 0, the system belongs to the dynamical class of model I indicating a dynamic 
phase transition at /? = 0. On the other hand, the persistence probabilities in a system of L spins 
saturate at a value P sa t{P, Li) = (/3 / L) a f '(/3) , where a remains constant for all ft 7^ supporting the 
existence of the dynamic phase transition at ft — 0. The scaling function /(/3) shows a crossover 
behaviour with f(/3) = constant for /? << 1 and f(/3) oc /3~ a for f3 » 1. 
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The effect of noise on equilibrium behavior is well 
known, e.g., there are order-disorder phase transitions in- 
duced by thermal noise observed in many systems. Noise 
induced phase transitions may occur in dynamical sys- 
tems as well when the noise can drive the system from 
one dynamical class to another. These dynamical classes 
are often characterised by different dynamical exponents. 
In this paper, we study a case of such a dynamical phase 
transition in a very simple Ising like spin system. 

A dynamical model of Ising spins has been recently 
proposed in [l[ (which we refer to as model I henceforth) 
where the state of the spins may change in two situations: 
first when its two neighbouring domains have opposite 
polarity, and in this case the spin orients itself along the 
spins of the neighbouring domain with the larger size. 
This case may arise only when the spin is at the bound- 
ary of the two domains. A spin is also flipped when it 
is sandwiched between two domains of spins with same 
sign. Except for the rare event when the two neighbour- 
ing domains of opposite spins are of the same size, the 
dynamics in the above model is deterministic. This dy- 
namics leads to a homogeneous state of either all spin 
up or all spin down. Such evolution to absorbing homo- 
geneous states are known to occur in systems belonging 
to directed percolation (DP) processes, zero temperature 
Ising model, voter model etc 0,0]. 

Model I was introduced in the context of a social sys- 
tem where the binary opinions of individuals are repre- 
sented by up and down spin states. In opinion dynamics 
models, such representation of opinions by Ising or Potts 
spins is quite common Q. The key feature is the inter- 
action of the individuals which may lead to phase tran- 
sitions between a homogeneous state to a heterogeneous 
state in many cases 

Model I showed the existence of novel dynamical be- 
haviour in a coersening process when compared to the 
dynamical behaviour of DP processes, voter model, Ising 
models etc In this work, we have introduced 

stochasticity in the dynamics of Model I to see how it 
affects the coarsening process. 

Let d up and ddown be the sizes of the two neighbour- 
ing domains of type up and down of a spin at the do- 



main boundary (excluding itself). In model I, probabil- 
ity P(up) that the said spin is up is 1 if d up > ddowni 
0.5 if d up — ddown and zero otherwise. In the sim- 
plest possible way to introduce stochasticity, one may 
take the probability of a boundary spin to be up as 
P(up) = d up /(d U p + ddown)- However, there is no pa- 
rameter controlling the stochasticity here and moreover, 
we find that the results are identical to the original model 
I. 

In order to introduce a noise like parameter which can 
be tuned, we next propose that the probability that a 
spin at the domain boundary is up is given by 

P{up) ae/K**.-**"."), (i) 

and it is down with probability 

P(down) cx e^*""""-**). (2) 

The normalised probabilities are therefore P(up) — 
exp/3A/(exp(/3A) + exp(-/3A)) and P(down) = 1 - 
P{up), where A = (d up - d down ). 

Obviously, (3 — > 00 cooresponds to model I while let- 
ting (3 = we have equal probabilities of the up and down 
states, making it equivalent to the zero temperature dy- 
namics of the nearest neighbour Ising model. Since the 
equilibrium states for the extreme values /3 — > 00 and 
(3 = are homogeneous (all up or all down states), it is 
expected that for all values of f3 they will be remain so 
as is indeed the case. 

As far as dynamics is concerned, we investigate primar- 
ily the time dependent behaviour of the order parameter 
and the persistence probability. In the one dimensional 
chain of length L, the order parameter is the conven- 
tional magnetisation given by M — w here 
L up (Ldown) is the number of up (down) spins in the 
system and L = L up + Ldowm the total number of spins. 
The average fraction of domain walls D w , which is the 
average number of domain walls divided by L is also stud- 
ied. D w is identical to the inverse of average domain size. 
The dynamical evolution of the order parameter and frac- 
taion of domain walls is expected to be governed by the 
dynamical exponent z; M ex. £ _1 /( 2z ) and D w ~ t~ z [llj ]. 
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The persistence probability of a spin is the probability 
that it remains in its original state upto time t [10|. It 
has been shown to have a power law decay in many sys- 
tems with an associated exponent 9. To obtain both the 
exponents 9 and z in finite systems of dimension L from 
the persistence probability, the following scaling form is 
often used [12] 

P{t 1 L)=t- e f{L/t 1 l z ). (3) 

Another exponent, a = 6z, is associated with the satura- 
tion value of the persistence probability at t —> oo when 
P sat (L) = P{t^ oo, L) cx L~ a [3. 

In model I, it was numerically obtained that 9 ~ 0.235 
and z ~ 1.0 giving a ~ 0.235, while in the one dimen- 
sional Ising model 9 — 0.375 and z — 2.0 (exact results) 
giving a — 0.75. It is clearly indicated that model I and 
the Ising model belong to two different dynamical classes. 
By introducing the parameter p one can therefore expect 
a transition from the Ising to the model I dynamical be- 
haviour at some specific value of p. 

With respect to model I, f3 — is the maximum noise 
and its inverse may be thought of an effective tempera- 
ture. On the other hand, from the Ising model viewpoint, 
P plays the role of noise. However it is not equivalent to 
thermal fluctuations which can affect the state of any 
spin. With P, flipping of spins can still occur at the do- 
main boundaries only. Hence, even with this noise, the 
equilibrium behaviour is not disturbed for any value of P 
(even for ft — > oo which corresponds to model I) while in 
contrast, any non-zero temperature can destroy the order 
of a one dimensional Ising model. 

It is useful to show the snapshots of the evolution of the 
system over time for different P (Fig.l): to be noted is 
the fact that for any non-zero P, the system equilibriates 
very fast compared to the Ising limit P = 0. 

In the simulations, we have generated systems of size 
1000 < L < 10000 with a mininum of 1000 initial con- 
figurations for the maximum size in general. Only for 
P = 0, the Ising limit, in which case the time taken to 
reach equilibrium is order of magnitude higher than that 
for any nonzero P, smaller systems have been simulated 
in some calculations. Depending on the system size and 
time to equilibriate, maximum iteration times have been 
set. Random sequential updating process has been used 
to control the spin flips. 

On introducing /3, we notice that well away from the 
Ising limit p = 0, the dynamics gives z ~ 1.0 and 
9 ~ 0.235 as in model I. However, as P is made less than 
O(10 _1 ), the behaviour of the relevant dynamic quan- 
tities deviate from a simple power law behaviour. For 
example, the magnetisation shows an initial slow varia- 
tion with time followed by a rapid growth before reaching 
saturation for values of P < 0.1 (Fig. 2). It is difficult 
to fit a power law in either regime. This is true for the 
domain wall fraction decay as well (not shown). In fact, 
the rapid growth of magnetisation at later times is ap- 
parently even faster than i 1 / 2 , that obtained for model 
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FIG. 1: (Color online) Snapshots of the system in time for 
different values of @ = 0.0,0.005,0.1 and /3 — > oo (top to 
bottom) showing that for any non-zero P, the system equilib- 
riates towards a homogeneous state much faster compared to 
the P = case. These snapshots are for a, L — 100 system. 



I (e.g., for p = 0.001). From the snapshots of the sys- 
tem for P very close to zero, it is seen that for the first 
few steps the system has a behaviour similar to the Ising 
model (P = 0). This explains the slow growth of magneti- 
sation initially. However, as soon as a domain shrinks in 
size compared to its adjacent one, any non-zero P makes 
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FIG. 2: (Color online) Magnetisation as a function of time is 
shown for /3 = 0, 0.001 and 0.01. The two starightlines in the 
log log plot have slopes corresponding to model I (0.5) and 
Ising model (0.25). The /3 = result is for L = 2000 while 
the others are for L — 10000. 



FIG. 3: (Color online) The values of the time to equilibriate 
t eq scaled by the sytem size shows that for large /3 there is 
a nice collapse. For small /3, there are deviations from the 
collapse which decrease with the system size. Inset shows 
that t eq scales as L z for the limiting values /3 = with z — 2 
and for f3 — > oo, 8 = 1. 



it vanish very rapidly. However, it will be wrong to in- 
fer that the coersening process takes place faster than in 
model I, because in comparison, in model I, the system 
cquilibriates in times much lesser than that for any finite 
P (see Fig. 1). 

The question remains therefore whether and how one 
can obtain an estimate of z for /3 — > 0. Since a direct fit- 
ting fails, we try an indirect method. The average time 
t eq to reach the equilibrium state can be estimated from 
the time the magnetisation reaches a value unity. t eq is 
shown to scale as L z in Ising model with z — 2 and for 
P — » oo, t eq scales as L z with z = 1 (inset of Fig. 3). 
Hence we plot t eq /L against P for different L and find 
an interesting result. For values of j3 greater than 0.01, 
it shows a nice collapse, indicating z = 1 here. As j3 
decreases, the deviation from a collapse starts appear- 
ing, it getting more pronounced for smaller values of /?. 
However, at the same time, we notice that the deviation 
from a scaling t eq oc L decreases for larger values of L 
suggesting that the collapse as j3 — > will improve with 
the system size. Thus we conclude that the exponent z 
equals unity in the thermodynamic limit for any non-zero 
value of p. The deviations from the scaling as j3 — > is 
simply a finite size effect. 

Hence from the above behaviour we conclude that the 
model I behaviour is valid for any finite j3 and a dynamic 
transition takes place exactly at = 0. 

Next we focus on the persistence data. Once again, 
as fi — >• 0, it is difficult to fit a unique power law to the 
persistence probability (Fig. 4). Here, in consistency 
with the magnetisation results, we find an initial decay 
of persistence quite fast and a late variation compara- 
tively slower. The initial variation can be fitted to a 
power law and an estimate of 9 made this way shows a 
tendency to continuously vary towards the P = value, 
i.e., 0.375. However, 9 is not to be obtained from the 
early time behaviour and there is definitely a crossover 
to a different behaviour in later times before the persis- 
tence reaches saturation. Therefore determining 9 from 
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FIG. 4: (Color online) The persistence probability against 
time is shown for different values of /3 (/? = 100, 0.1, 0.01 and 
0.005 (from top to bottom); for small j3 the slope cannot be 
uniquely determined. The two straightlines in the log log plot 
have slopes corresponding to the model I (0.235) and Ising 
model (0.375). 



the initial variation is not a correct approach. 

We even try to obtain a collapse by plotting P(t)/t~ e 
against L/t z using trial values of z and 9 as in [l|, [l3[ , 
but for P — » 0, no collapse for large L/t z , i.e., for small t, 
can be obtained, confirming once again that the determi- 
nation of 9 is not possible in a straightforward manner. 

We next try to find out whether the scaling law 
Psat(L) oc L~ a is valid for finite values of p. When we 
plot the saturation values of persistence against differ- 
ent system sizes, we do find nice power law fittings and 
hence estimates of a can be made (Fig. 5). We find 
that a varies between 0.22 and 0.23 with no systematics 
indicating that it is independent of p. This once again 
supports the fact that there is a transition at P — as a 
has a known value (0.75) much larger for P = 0. 
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FIG. 5: (Color online) The saturation values of the persis- 
tene probability shows a variation L~ a for values of /? = 
10.0, 0.1, 0.01 and 0.001 (from top to bottom) with a = 0.225. 
Inset shows that the scaled saturation values P 3a t / (P / L) a — 
/(/?) varies as fi~ a for large ft. 



Although a shows no dependence on /3, the saturation 
values of the persistence probability show an interesting 
dependence on /3: for j3 >> 1, it is independent of (3 while 
for small values of (3 it has a power law variation. We 
in fact find that the scaled variable P sa t/(^) a with a = 
0.225 shows a collapse when plotted against f3 suggesting 
a scaling form: 



P sat {L,p) = {p/L) a f{p). 



(4) 



The fact that for large values of (3, the saturation values 
are independent of (3 suggests that f(/3) varies as (3~ a 



here. We indeed find this kind of a behaviour with /(/?) = 
constant for (3 « 1 and f(J3) oc (3~ a for (3 » 1 (see 
inset of Fig. 5). 

We thus find that the effect of the noise parameter f3 
is to cause a dynamic phase transition at /3 = showing 
that the behaviour of model I is indeed very robust. On 
the other hand, with respect to the Ising model, although 
the effect of noise is not comparable to thermal fluctua- 
tions as far as order-disorder transitions are concerned, 
it does induce a dynamic phase transition at p = 0. The 
signature of the dynamical phase transition is seen in the 
variation of the dynamical quantities as the j3 = point 
is approached, there are also strong finite size effects. 

One may raise the question as to what happens if /3 
is made negative. As expected, the system goes to a 
disordered state for any non-zero /? accompanied with 
exponential decay of persistence probability. In the spin 
picture, a negative value of (3 does not correspond to any 
physical model, but in terms of domain wall movement, 
one has a system of mutually repulsive random walkers 
when /3 < 0. The random walkers tend to move away 
from their nearest neighbours and therefore cannot an- 
nihilate each other but remain mobile all the time de- 
stroying the persistence of the spins. Such situations was 



seen to arise in spin systems like the ANNNI model [14 1 
also. The dynamic behaviour is therefore different for 
(3 < 0,(3 — and (3 > 0. So, allowing negative values of 
(3, one may say that there is a dynamic phase transition 
occurring at /3 = separating three different dynamical 
phases. 
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